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ABSTRACT 



We present the final results of a ROSAT PSPC program to study the soft X-ray 
emission properties of a complete sample of of low z quasars. This sample includes all 
23 quasars from the Bright Quasar Survey with z < 0.400, and N§^< 1.9 x 10 20 cm -2 . 
Pointed ROSATPSPC observations were made for all quasars, yielding high S/N spectra 
for most objects which allowed an accurate determination of the spectral shape. The 
following main results were obtained: 

1. The spectra of 22 of the 23 quasars are consistent, to within ~ 30%, with a single 
power-law model at rest-frame 0.2 — 2 keV. There is no evidence for significant 
soft excess emission with respect to the best fit power-law. We place a limit (95% 
confidence) of ~ 5 x 10 19 cm -2 on the amount of excess foreground absorption by 
cold gas for most of our quasars. The limits are ~lx 10 19 cm -2 in the two highest 
S/N spectra. 

2. The mean 0.2 — 2 keV continuum of quasars agrees remarkably well with an ex- 
trapolation of the mean 1050A-350A continuum recently determined by Zheng et 
al. (1996) for z > 0.33 quasars. This suggests that there is no steep soft component 
below 0.2 keV. 
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3. Significant X-ray absorption (r > 0.3) by partially ionized gas ("warm absorber") 
in quasars is rather rare, occurring for <^ 5% of the population, which is in sharp 
contrast to lower luminosity Active Galactic Nuclei (AGNs), where significant ab- 
sorption probably occurs for ~ 50% of the population. 

4. Extensive correlation analysis of the X-ray continuum emission parameters with 
optical emission line parameters indicates that the strongest correlation is between 
the spectral slope a x , and the H/3 FWHM. A possible explanation for this remarkably 
strong correlation is a dependence of a x on L/L Edd , as seen in Galactic black hole 
candidates. 

5. The strong correlations between a X) and L[ Q m ], Fe II/H/3, and the peak [O III] 
to H,3 flux ratio are verified. The physical origin of these correlations is still not 
understood. 

6. There appears to be a distinct class of "X-ray weak" quasars, which form ~ 10% 
of the population (3 out of 23), where the X-ray emission is smaller by a factor 
of 10-30 than expected based on their luminosity at other bands, and on their H/3 
luminosity. These may be quasars where the direct X-ray source is obscured, and 
only scattered X-rays are observed. 

7. Thin accretion disk models cannot reproduce the observed 0.2-2 keV spectral shape, 
and they also cannot reproduce the tight correlation between the optical and soft 
X-ray emission. An as yet unknown physical mechanism must be maintaining a 
strong correlation between the optical and soft X-ray emission. 

8. The H I/He I ratio in the high Galactic latitude ISM must be within 20%, and 
possibly within 5%, of the total H/He ratio of 10, which indicates that He in the 
diffuse H II gas component of the interstellar medium is mostly ionized to He II or 
He III. 

We finally note the intriguing possibility that although (a x ) in radio-loud quasars 
(— 1.15±0.14) is significantly flatter than in radio-quiet quasars (— 1.72±0.09) the X-ray 
emission may not be related to the presence of radio emission. The difference in (a x ) 
may result from the strong a x vs. H/3 FWHM correlation and the tendency of radio-loud 
quasars to have broader H/3. 

Subject headings: galaxies: active — quasars: emission lines — X-rays: galaxies — ISM: 
abundances 
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1. INTRODUCTION 

Quasars emit most of their power in the UV to soft 
X-ray regime. The PSPC detector aboard ROSAT al- 
lowed a significantly improved study of the soft X-ray 
emission of quasars compared with earlier missions 
(some of which were not sensitive below 2 keV) , such 
as HEAO-1, EINSTEIN, EXOSAT, and GINGA (e.g. 
Mushotzky 1984; Wilkes & Elvis 1987; Canizares & 
White 1989; Comastri et al. 1992; Lawson et al. 
1992; Williams et al. 1992; and a recent review by 
Mushotzky Done, & Pounds 1993). These earlier 
studies indicated that the X-ray emission above 1- 
2 keV is well described by a power law with a spectral 
slope a x — dlnf v /dlni/ of about —0.5 for radio loud 
quasars and about —1.0 for radio quiet quasars. Large 
heterogeneous samples of AGNs were recently studied 
using the ROSATPSPC by Walter & Fink (1993) and 
by Wang, Brinkmann, & Bergeron (1996). However, 
the objects studied in the papers mentioned above do 
not form a complete sample, and the available results 
may be biased by various selection effects which were 
not well defined a priori. In particular, most of the 
studied objects are nearby, intrinsically X-ray bright 
AGNs. 

To overcome the potential biases in existing studies 
we initiated a ROSATPSPC program to make an ac- 
curate determination of the soft X-ray properties of a 
well defined and complete sample of quasars, selected 
independent of their X-ray properties. This program 
was designed to address the following questions: 

1. What are the soft X-ray spectral properties of 
the low redshift quasar population? 

2. Are simple thin accretion disk models (e.g. Laor 
1990) able to fit the observed optical/UV/soft 
X-ray continuum? are other modifying mecha- 
nisms, such as a hot corona (e.g. Czerny & Elvis 
1987) required? Arc models invoking optically 
thin free-free emission possible (e.g. Ferland, 
Korista & Peterson, 1990; Barvainis 1993)? 

3. Do the observed soft X-ray properties display 
any significant correlations with other proper- 
ties of these quasars? Arc these correlations 
compatible with various models for the contin- 
uum and line emission mechanisms? 

Our sample includes all 23 quasars from the BQS 
sample (Schmidt & Green 1983) with z < 0.400, and 



N§f< 1.9 x 10 20 cm" 2 , where N§°} is the H I Galac- 
tic column density as measured at 21 cm. These selec- 
tion criteria allow optimal study of soft X-ray emis- 
sion at the lowest possible energy. The additional 
advantages of the BQS sample are that it has been 
extensively explored at other wavelengths (see Pa- 
per I for further details), and that it includes only 
bright quasars, thus allowing high S/N X-ray spectra 
for most objects. The sample selection criteria are 
independent of the quasar's X-ray properties, and we 
thus expect our sample to be representative of the 
low-redshift, optically-selected quasar population. 

Preliminary results from the analysis of the first 10 
quasars available to us were described by Laor et al. 
(1994, hereafter Paper I). Here we report the analysis 
of the complete sample which allows us to address the 
three questions posed above. The outline of the paper 
is as follows. In §2 we describe the observations and 
the analysis of the spectra. §3 describes the analy- 
sis of correlations between the soft X-ray properties 
and other continuum and emission line properties. In 
§4 we compare our results with other soft X-ray ob- 
servations and discuss some of the implications. We 
conclude in §5 with answers to the questions raised 
above, and with some new questions to be addressed 
in future studies. 

2. THE OBSERVATIONS AND ANALYSIS 
OF THE SPECTRA 

The complete sample of 23 quasars is listed in 
Table 1 together with their redshifts, ms and Mb 
magnitudes (calculated for H = 50 km s _1 , q=0.5), 
R, the radio to optical flux ratio, and N^j. The 
redshifts, ms, and Mb magnitudes are taken from 
Schmidt & Green (1983), R is taken from Kcllcrmann 
et al. (1989). Note that 4 of the 23 quasars in our 
sample are radio loud (defined here as R > 10). 

The Galaxy becomes optically thick below 0.2 keV 
for the typical high Galactic column of 3 x 10 20 cm -2 
(Dickey & Lockman 1990; Morrison & McCammon 
1983), and accurate values of A^y 1 are therefore cru- 
cial even for our low A^fy 1 sample. The A^fy 1 values 
given in column 7 of Table 1 are taken from Elvis, 
Lockman & Wilkes (1989), Savage et al. (1993), 
Lockman & Savage (1995), and the recent extensive 
measurements by Murphy et al (1996). All these 
measurements of iVjfy were made with the 140 foot 
telescope of the NRAO at Green Bank, WV, us- 
ing the "bootstrapping" stray radiation correction 
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method described by Lockman, Jahoda, & McCam- 
mon (1986), which provides an angular resolution of 
21', and an uncertainty of AA^fy=l x 10 19 cm -2 (and 
possibly lower for our low A^h quasars). This uncer- 
tainty introduces a flux error of 10% at 0.2 keV, 30% 
at 0.15 keV, and nearly a factor of 2 at 0.1 keV. Thus, 
reasonably accurate fluxes can be obtained down to 
- 0.15 keV. Note that Murphy et al. (1996) includes 
accurate measurements for about 220 AGNs, 

including most AGNs observed by ROSAT, which 
would be very useful for eliminating the significant 
systematic uncertainty in the PSPC spectral slope 
which must be present when a low accuracy A^ 1 is 
used. 

Table 2 lists the PSPC observations of all the 
quasars. For the sake of completeness we include 
also the 10 quasars already reported in Paper I. All 
sources were detected, and their net source counts 
range from 93 to 38,015, with a median value of about 
1900 counts. The PROS software package was used 
to extract the source counts. Table 2 includes the 
exposure times, the dates of the observations, the 
net number of counts and their statistical error, the 
count rate, the radius of the circular aperture used 
to extract the source counts, the offset of the X-ray 
position from the center of the PSPC field of view, 
the ROSAT sequence identification number, and the 
SASS version used for the calibration of the data. All 
objects, except one, are typically within 15" of the 
center of the PSPC field of view, so all the identifica- 
tions are secure. The one exception is PG 1440+356 
where the pointing was offset by 40' from the posi- 
tion of the quasar (Gondhalekar et al. 1994). Note 
that the exposure times are uncertain by about 4% 
due to a number of possible systematic errors, as de- 
scribed by Fiore et al. (1994). The typically large 
number of counts for each object allows an accurate 
determination of the spectral slope for most objects, 
as described below. 

Model fits to the extracted number of source counts 
per pulse invariant (PI) channel, N°^, were carried 
out using the XSPEC software package. PI chan- 
nels 1-12 of the original 256 channels spectra (E < 
0.11 keV), were ignored since they are not well cali- 
brated and are inherently uncertain due to the large 
Galactic optical depth. The January 1993 PSPC cal- 
ibration matrix was used for observations made after 
1991 Oct. 14, and earlier observations were fit with 
the March 1992 calibration matrix. The best fit model 
parameters are obtained by \ 2 minimization. Nearby 



channels were merged when < 10. A 1% error 
was added in quadrature to the statistical error in 
, to take into account possible systematic calibra- 
tion errors (see Paper I for more complete details). 

2.1. A Single Power-Law 

As in Paper I, we fit each spectrum with a single 
power-law of the form Je = e~ NH - aE foE a * , where /e 
is the flux density oe is the absorption cross section 
per H atom (Morrison & McCammon 1983), /o is the 
flux density at 1 keV, and E is in units of keV. We 
make three different fits for each object, with: 1. A^h 
a free parameter, 2. N U =N^. 3. N n =N§f, and 
0.47 < E < 2.5 keV, i.e. using only the hard ROSAT 
band (channels 12 — 34 of the rebinned 34 channels 
spectra). A comparison of fits 1 and 2 allows us to 
determine whether there is evidence for a significant 
intrinsic absorption excess or emission excess relative 
to a single power-law fit with Nu=N^ l f. This com- 
parison also allows us, as further shown in §5.3.3, to 
determine whether the 21 cm measurement of A^f^ 1 is 
a reliable measure of the Galactic soft X-ray opacity. 
A comparison of fits 2 and 3 allows us to look for a 
dependence of the power-law slope on energy. 

Table 3 provides the results of fits 1-3 described 
above. The table includes the 13 objects not reported 
in Paper I, and 6 objects from paper I for which we 
only now have the accurate Afjfy values. For each 
fit we give the best fitting spectral slope a x , the nor- 
malization of the power-law flux at 1 keV (/o), the 
best fitting Ah, the \ 2 01 the fit (x| t ), the number 
of degrees of freedom (dof), and the probability for 
X 2 > Xfit- The errors Aa x and AA^h in fit 1 were 
calculated by making a grid search for models with 
Ax 2 — 2.30, as appropriate for 1 a confidence level 
for two interesting parameters (e.g. Press et al. 1989). 
The error on the slope Aa. x in fits 2 and 3 is calculated 
by requiring A\ 2 = 1.0 (i.e. 68% for one interesting 
parameter). We neglect the effect of AATjfy on Aa x 
in fits 2 and 3 since we use accurate for all ob- 

jects. 

The observed and best fit spectra for the 13 quasars 
not reported in paper I are displayed in Figure 1. 
There are 3 panels for each object. The upper panel 
displays the observed count rate per keV as a func- 
tion of channel energy, the histogram represents the 
expected rate from the best fit power-law model, N°£, 
with A^h a free parameter (fit 1). The middle panel 
displays A/a, where A = N°£ - N^ od , and a is the 
standard error in N°£. This plot helps indicate what 
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features in the spectrum are significant. The lower 
panel displays the fractional deviations from the ex- 
pected flux, or equivalcntly A/N™ d , which indicates 
the fractional amplitude of the observed features. 

As shown in Table 3, in all 13 objects the sim- 
ple power-law model with (fit 2) provides an 
acceptable fit (i.e. Pr(x 2 > xlt) > 0.01). Note in 
particular the spectrum of PG 1116+215, which de- 
spite the very high S/N available (24,272 net counts), 
shows no deviations from a simple power law above a 
level of <; 10%. In Paper I a simple power-law model 
could not provide an acceptable fit to three of the 10 
quasars, though in two of them the apparent features 
could not be fit with a simple physical model, and in 
one of them this may be due to calibration errors (see 
§4.1). 

As mentioned above, a comparison of the free TVh 
fit (fit 1) with the N§f fit (fit 2) allows us to look for 
evidence for an absorption or an emission excess. We 
measure the statistical significance of the reduction in 
Xjj t with the addition of Nn as a free parameter using 
the F test (Bevington 1969). In PG 1440+356 we find 
a significant reduction with Pr= 7.5 x 10~ 4 (F=12.98 
for 26 dof), where Pr is the probability that the re- 
duction in x 2 is not statistically significant (calculated 
using the FTEST routine in Press et al. 1989). The 
./Vh obtained in fit 1 suggests intrinsic absorption of 
about 5 x 10 19 cm~ 2 above the Galactic absorption. 
However, unlike all other objects, PG 1440+356 was 
observed significantly off axis (see Table 2), and some 
small systematic calibration errors may be present 
there. Note also that the xlt °f the fit with -/V^ff 
(Pr=0.05) is still acceptable. We therefore cannot 
conclude that an extra absorber must be present in 
PG 1440+356. Marshall et al. (1996) found a very 
steep slope (a = -4.7+0.65) in PG 1440+35 at 0.1- 
0.15 keV (80-120A) using the Extreme UV Explorer. 
There is no indication for such a component in the 
PSPC spectrum below 0.15 keV (Fig.lc). However, 
given the very low sensitivity of the PSPC below 
0.15 keV, such a steep soft component may still be 
consistent with the PSPC spectrum. In the other 12 
objects the free TVh fit does not provide a significant 
improvement (i.e. Pr> 0.01), and thus there is no 
clear evidence for either intrinsic absorption, or low 
energy excess emission above a simple power-law. 

Figure 2 compares the Galactic Nn deduced from 
the accurate 21 cm measurements with the best fit 
X-ray column deduced using the free Nh fit. The 
straight line represents equal columns. The x 2 of 



the 7Vh(21 cm) =W H (X-ray) model is 31.9 for 22 dof 
(PG 1001+054 was not included because of the low 
S/N), which is acceptable at the 8% level. This re- 
sult demonstrates that there is no significant excess 
absorption over the Galactic value in any of our ob- 
jects. It is interesting to note that in our highest S/N 
spectra, PG 1116+215 and PG 1226+023, iV H (X-ray) 
is determined to a level of 0.8 — 1 x 10 19 cm 2 , and is 
still consistent with ATh(21 cm), indicating that both 
methods agree to better than 10%. 

The average hard ROSAT band (0.5-2 keV) slope 
for the complete sample is —1.59 ± 0.08 (excluding 
PG 1114+445 which is affected by a warm absorber, 
and PG 1001+054 and PG 1425+267 where the S/N 
is very low). This slope is not significantly differ- 
ent from the average slope for the full ROSAT band, 
-1.63 + 0.07. 

Spectral fits to the PSPC data of some of the 
objects in our complete sample have already been 
reported by Gondhalekar et al. (1994), Ulrich & 
Molcndi (1996), Rachcn, Mannheim, & Bicrmann 
(1996), and Wang, Brinkmann, & Bergeron (1996). 
The results of the single power-law fit with a free N-r 
in these papers are all consistent with our results for 
the overlapping objects. The only discrepancy is with 
PG 1444+407 where although both us and Wang et 
al. find a similar slope, Wang et al. find evidence 
for absorption while we find no such evidence. For 
a simple power-law fit to PG 1444+407 Wang et al. 
find x 2 = 26 for 20 dof which is acceptable only at 
the 17% level (~ 1.4er level), while we find for such 
a fit x 2 = 14-7 for 20 dof, which is acceptable at the 
80% level. 

As discussed in paper I (§5.1.3), the difference in 
spectral slopes at hard (2-10 keV) and soft X-rays 
raises the possibility that a x may be changing within 
the PSPC band itself. The individual spectra are well 
fit by a simple power law, and thus any spectral cur- 
vature must be consistent with zero. Stronger con- 
straints on the spectral curvature may be obtained by 
measuring the average curvature parameter (/3, de- 
fined in Paper I) for the complete sample since the 
random error in the mean is smaller by \/N <~ 5 
than the random error for individual objects. Un- 
fortunately, the PSPC calibration uncertainty at low 
energy, discussed in Paper I, introduces a systematic 
error in [3 (which obviously does not cancel out as 
VN), and as shown in paper I, does not allow a re- 
liable determination of the curvature parameter. We 
therefore did not try to constrain the spectral curva- 
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ture parameter in this paper. 

3. CORRELATION ANALYSIS 

Table 4 presents 8 of the 12 rest-frame contin- 
uum parameters and 7 of the 18 emission line pa- 
rameters used for the correlation analysis. The spec- 
tral slopes arc defined as (flux subscript indicates 
logi/(Hz)): a Q = log(/i 5 //i 4 .5)/0.5 (1-0.3 /zm), a ox = 
log(/i 7 .68 5 //i 5 )/2.685 (3000 A-2 keV), and a os = 
log(/i 6 .86i//i 5 )/1.861 (3000 A-0.3 keV). The X-ray 
continuum parameters are from fit 2, and from Pa- 
per I. The near IR and optical continuum parame- 
ters are taken from Neugebauer et al. (1987). The 
emission line parameters were taken from Boroson & 
Green (1992). Luminosities were calculated assuming 
if =50 km s _1 and q = 0.5. 

The four additional continuum parameters not 
presented in Table 5 for the sake of brevity are 
«irx = log(/i 6 .86i//i4.2 5 )/2.611 (1.69 ^m-2 keV), 
"irs = log(/i 7 .68 5 //i4.25)/3.435 (1.69 M m-0.3 keV), ra- 
dio luminosity (Kcllcrmann et al. 1989), and 1 /im lu- 
minosity (Neugebauer et al. 1987). The 11 emission 
line parameters not detailed in Table 5 are: [O III] 
EW, Fe II EW, H/3 EW, He II EW, [O III] /FT? and 
He II/H/3 flux ratio, [O III] peak flux to H/3 peak 
flux ratio, radio to optical flux ratio, and the H/3 as- 
symetry, shape, and shift parameters. All these 11 
parameters are listed in Table 2 of Boroson & Green 
(1992). 

The significance of the correlations was tested us- 
ing the Spearman rank-order correlation coefficient 
(rs) which is sensitive to any monotonic relation be- 
tween the two variables. A summary of the main cor- 
relation coefficients and their two sided significance is 
given in Table 5. 

3.1. The Significance Level 

In Paper I the correlation analysis was carried out 
using 10 objects, and only relatively strong correla- 
tions (rs > 0.76) could be detected at the required 
significance level (Pr< 0.01). Here, with 23 objects, a 
Pr< 0.01 corresponds to rs > 0.52, and we can thus 
test for the presence of weaker correlations, and check 
whether the correlations suggested in Paper I remain 
significant. We have searched for correlations among 
the 12 continuum emission parameters, and between 
these 12 parameters and the 18 emission line param- 
eters listed above, which gives a total of 294 different 
correlations. One thus expects about 1 spurious cor- 



relation with Pr< 3.4 x 10 3 in our analysis, and for 
a significance level of 1% one would now have to go 
to Pr< 3.4 x 10~ 5 , rather than Pr< 1 x 10~ 2 . How- 
ever, we find that there are actually 42 correlations 
with Pr< 3.4 x 10~ 3 , rather than just one, in our 
sample. Thus, the probability that any one of them 
is the spurious one is only 2.4%, and the significance 
level of these correlations is reduced by a factor of 7 
(= 0.024/0.0034), rather than a factor of 300. Below 
we assume that correlations with Pr< 1 x 10~ 3 are 
significant at the 1% level (there are 30 correlations 
with Pr< 1 x 10~ 3 , versus an expected number of 
0.3). Thus, given the large number of correlations we 
looked at, we can only test reliably for correlations 
with rs > 0.64 (which corresponds to Pr< 1 x 10~ 3 
for 23 data points). 

3.2. The Near IR to X-Ray Energy Distribu- 
tion 

A comparison of the rest-frame spectral energy dis- 
tributions of all 23 quasars is shown in Figure 3. The 
3-0.3 /im continuum is from Neugebauer et al. (1987), 
and the 0.2-2 keV continuum is from paper I and 
from this paper. The upper panel shows the abso- 
lute luminosities, and the lower 2 panels the lumi- 
nosity normalized to unity at logf = 14.25 for radio 
quiet quasars and for radio loud quasars. Note the 
relatively small dispersion in the normalized 0.3 keV 
(logzv = 16.861) luminosity. The outlying objects are 
labeled. PG 1626+554 is the only object where a 
steep a x is clearly associated with a strong soft ex- 
cess (relative to the near IR flux). In other objects 
a steep a x tends to be associated with a low 2 keV 
flux. This trend is also suggested by the presence of 
a marginally significant correlation between a x and 
otox (rs — 0.533, Pr= 0.0089, see below), and the ab- 
sence of a significant correlation between a x and a os 
(r s = -0.201, Pr= 0.36). 

The X-ray luminosity distribution appears to be bi- 
modal with two quasars, PG 1001+054 and PG 1411+ 
442, being a factor of 30 weaker than the mean radio 
quiet quasar. These two quasars appear to form a 
distinct group of 'X-ray weak quasars'. The statistics 
for the radio loud quasars (RLQ) arc much poorer, 
and there is no well defined mean, but PG 1425+267 
may be a similar X-ray weak RLQ. 
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3.3. Correlations with Emission Line Proper- 
ties 

Figure 4 presents the correlations between the hard 
X-ray luminosity, L 2 k e v (logf = 17.685), or the soft 
X-ray luminosity, io.3 kcV> with the luminosities of 
Up, [O III], He II, or Fc II. The value of r s and the two 
sided significance level (Pr) of rs are indicated above 
each panel. Upper limits were not included in the 
correlations. Thus, the actual correlations for He II, 
where there are 5 upper limits, are likely to be smaller 
than found here (there is only one upper limit for 
[O III] and Fe II, and none for H/3). Excluding He II, 
the X-ray luminosity is most strongly correlated with 
(rs = 0.734, Pr= 6.6 x 10~ 5 ). We note in passing 
that Lr/3 has an even stronger correlation with the 
luminosity at 3000 k(r s = 0.866, Pr= 9 x 10~ 8 ) and 
with the near IR luminosity at 1 /im (rs = 0.810, 
Pr= 2.7 x 10~ 6 ). 

The position of the X-ray weak quasars is marked 
in Fig.4. Both PG 1001+054 and PG 1411+442 ap- 
pear to have an X-ray luminosity weaker by a factor 
of about 30 compared to other quasars with similar 
Lh/3- PG 1425+267 is also weaker by a factor of <~ 10 
compared with the other RLQ. These ratios are the 
same as those found above in §3.2, based on the spec- 
tral energy distribution. 

Figures 5a-d displays various emission parameters 
which correlate with a x (as obtained with Nh^N^). 
The FWHM of H/3, L [Q in], the Fe II/H/3 flux ratio, 
and the ratio of [O III] peak flux to H^ peak flux (as 
defined by Boroson & Green) arc the emission line 
parameters which correlate most strongly with a x . 
As found in Paper I, all the a x versus emission line 
correlations become significantly weaker when we use 
a x obtained with the free Ah fit. 

The X-ray weak quasars are labeled in the a ox vs. 
a x correlation in Figure 5e. As expected they have a 
steeper than expected a ox for their a x . The last pa- 
rameter shown in Figure 5f is I.HL\^ 2 25 Av~ 2 , where 
Av =H/3 FWHM. This parameter is related, under 
some assumptions, to the luminosity in Eddington 
units, as further discussed in §4.7. 

4. DISCUSSION 

4.1. The Soft X-Ray Spectral Shape 

We find an average spectral index (a x ) — —1.62 + 
0.09 for the complete sample of 23 quasars, where 
the error here and below is the uncertainty in the 



mean. This slope is consistent with the mean slope 
(&x) = —1-57+0.06 which we found for the subsample 
of 24 quasars out of the 58 AGNs analyzed by Walter 
& Fink (1993, the other 34 AGNs in their sample are 
Seyfert galaxies as defined by Veron-Cetty & Veron 

1991) . A similar average slope of —1.65 ± 0.07 was 
found by Schartel et al. (1996) for 72 quasars from 
the LBQS sample detected in the i?0,S^4Tall sky sur- 
vey (RASS). Puchnarewicz et al. (1996) find a sig- 
nificantly flatter mean slope, (a x ) — —1.07+ 0.06, in 
a large sample of 108 soft X-ray selected (0.5-2 keV) 
AGNs. Part of the difference is related to the exclu- 
sion of counts below 0.5 keV, which selects against 
steep a x quasars, but this bias cannot explain the 
much flatter a ox in their sample. As discussed by 
Puchnarewicz et al., their sample appears to include 
a large proportion of highly reddened quasars (see 
further discussion in §4.2). 

RLQ are known to have a flatter a x than radio 
quiet quasars (RQQ) at energies above the PSPC 
band (e.g. Wilkes & Elvis 1987, Lawson et al. 1992). 
We find (a x ) = -1.72 ± 0.09 for the 19 RQQ, and 
(a x ) = -1.15 + 0.14 for the 4 RLQ in our sample. We 
find a similar trend using the Walter & Fink quasar 
data, where (a x ) — —1.61 ± 0.08 for the RQQ and 
(a x ) = -1.36 ± 0.08 for the RLQ. A similar differ- 
ence between RQQ and RLQ was found by Ciliegi 
& Maccacaro (1996) in PSPC spectra of a sample 
of 63 AGNs extracted from the EINSTEIN extended 
medium sensitivity survey sample. We therefore con- 
clude that the trend observed at harder X-rays also 
extends down to the 0.2-2 keV band. Schartel et al. 
(1995) stacked PSPC images of 147 RQQ and 32 RLQ 
finding for the sum images (a x ) = — 1.65 ± 0.18 for 
RQQ and (a x ) = -1.00 + 0.28 for the RLQ. However, 
the mean redshift of their objects is ~ 1.3 and thus 
their results apply to the ~ 0.45 — 4.5 keV band. 

As discussed in Paper I, the ROSAT PSPC indi- 
cates a significantly different soft X-ray spectral shape 
for quasars compared with earlier results obtained by 
the EINSTEIN IPC and EXOSAT LE+ME detec- 
tors (e.g. Wilkes & Elvis, 1987; Masnou ct al. 1992; 
Comastri et al. 1992; Saxton et al. 1993; Turner & 
Pounds 1989; Kruper, Urry & Canizares 1990). In 
particular, earlier missions suggested that the hard 
X-ray slope (Lawson et al. 1992; Williams et al. 

1992) extends down to <~ 0.5 keV with a steep rise at 
lower energy. Here we find that the 0.2-2 keV spec- 
trum is fit well by a single power-law with Galactic 
absorption. This indicates that: 1). the break be- 
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tween the soft and hard X-ray slope must occur well 
above 0.5 keV, 2) the break must be gradual, and 
3) there is no steep soft component with significant 
flux down to <~ 0.2 keV. ASCA observations of two of 
the quasars in our sample, 3C 273 by Yaqoob et al. 
(1994), and PG 1116+215 by Nandra et al. (1996) 
find, as expected, significantly flatter spectra above 
2 keV. However, the exact break energy cannot be 
accurately determined from the ASCA spectra due to 
likely calibration uncertainties below 1 keV. 

As mentioned in Paper I, the different EINSTEIN 
IPC and EXOSAT LE+ME results may be traced 
back to the combined effect of the lower sensitivity 
of these instruments below ~ 0.5 keV, and possibly 
some calibration errors. Small systematic errors in 
the PSPC response function appear to be present be- 
low 0.2 keV (Fiore et al. 1994), and this instrument 
is thought to be significantly better calibrated at low 
energy than earlier instruments. 

No significant spectral features are present in the 
PSPC spectra of all 13 additional quasars reported 
here, indicating that intrinsic features must have an 
amplitude of less than 10-20%. Note, in particular 
the high S/N spectrum of PG 1116+215, where the 
number of counts is about 12 times the median sample 
counts, yet this spectrum is still consistent with a 
simple power-law. For the complete sample we find 
that only 1 quasar, PG 1114+445, has a significant 
physical feature which is well described by a warm 
absorber model. In two other quasars, PG 1226+023 
and PG 1512+370, there are significant features below 
0.5 keV (paper I). In the case of PG 1512+370 the 
features are at a level of <~ 30%, and in PG 1226+023 
they are at a level <;10% and may well be due to small 
calibration errors. This result is consistent with the 
result of Fiore et al. (1994) who found that a simple 
power-law provides an acceptable fit to the individual 
spectra of six high S/N PSPC quasar spectra. 

A composite optical to hard X-ray spectral energy 
distribution for RLQ and RQQ is displayed in Figure 
6. To construct it we used the mean £14.25 (Table 
5), the mean a , the mean a x , and the mean a ox in 
our sample. We excluded from the mean the three 
X-ray weak quasars, and PG 1114+445, where a x is 
highly uncertain due to the presence of a warm ab- 
sorber. The mean spectra were extended above 2 keV 
assuming a slope of —1 for RQQ and —0.7 for RLQ. 
The Mathews & Ferland (1987, hereafter MF) quasar 
energy distribution is also displayed for the purpose 
of comparison. The MF shape assumes a steep soft 



component with a break to the hard X-ray slope above 
0.3 keV, and it therefore significantly underestimates 
the soft X-ray flux at - 0.2 - 1 keV. 

RLQ tend to be somewhat stronger hard X-ray 
sources than RQQ. This trend, together with the flat- 
ter X-ray slope of RLQ was interpreted by Wilkes & 
Elvis (1987) as possible evidence for a two component 
model. In this interpretation RLQ have the same hard 
X-ray component with a x <~ —1, as in RQQ, with 
an additional contribution from a flatter a x <~ —0.5 
component, making their overall X-ray emission flat- 
ter and brighter. The additional X-ray component in 
RLQ could be related to the radio jet, e.g. through in- 
verse Compton scattering. The composite spectrum 
suggests that although RLQ are brighter at 2 keV, 
they may actually be fainter at lower energy because 
of their flatter a x . The RLQ composite is based only 
on four objects and is therefore rather uncertain. In 
addition, the results of Sanders et al. (1989 §111. c) 
suggest that RLQ in the PG sample are about twice 
as bright at 2 keV compared with RQQ of similar 
optical luminosity, rather than the <~ 30% found for 
the composite, thus the difference in PSPC a x would 
imply a smaller difference at 0.2 keV than shown in 
the composite. If RLQ are indeed weaker than RQQ 
at 0.2 keV then the two component model suggested 
above would not be valid, and RLQ need to have a 
different X-ray emission process, rather than just an 
additional component. 

The difference between RLQ and RQQ may actu- 
ally be unrelated to the radio emission properties, as 
discussed in §4.4. 

Figure 6 also displays a simple cutoff power-law 
model of the form L v cx v a ° e - hv / kT ™t w j t ^ Uo _ 
—0.3 and T cut — 5.4 x 10 5 K. This is an alternative 
way to interpolate between the UV and soft X-ray 
emission, and it is also a reasonable approximation 
for an optically thick thermal component. The lack 
of a very steep low energy component down to 0.2 keV 
allows us to set an upper limit on T cut . The upper 
limit is set using a and a os i the slope from 3000 A to 
rest frame 0.15(1 + z) keV (the lowest energy where 
the Galactic absorption correction error< 30%), given 

by 

cw = [2.685a o:E -(17.685-log!vK]/log(!v/10 15 ), 

where \ogVs' — 16.560 + log(l + z). 
The upper limit on the cutoff temperature T^ t is rc- 
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lated to the spectral slopes by 

, = 4.8 x 10- n (^ -10 15 )loge 

(a o -a os 01og(^/10 15 ) ' 

We find a rather small dispersion in T c u u ' t with (T c u u ' t ) = 
(5.5 ± 2.6) x 10 5 K, averaged over the complete sam- 
ple (Table 4), which corresponds to a cutoff energy 
of 47 cV, or about 3.5 Ryd. This value of T^ t cor- 
responds very closely to the far UV continuum shape 
assumed by MF (see Fig. 6). 

Walter et al. (1994) fit such a cutoff model di- 
rectly to six quasars and Seyfert galaxies finding 
(S cut ) = 63 ± f 2 eV, or (T cut ) = (7.3 ± 1.4) x 10 5 K, 
while Rachen, Mannheim & Biermann (1996) find us- 
ing such a model (T cut ) = (6.3 ± 2.3) x 10 5 K for 7 
quasars and Seyfert galaxies. These values are consis- 
tent with our results. The small dispersion in T cut re- 
flects the small dispersion in a os in our sample, which 
is in marked contradiction with the dispersion pre- 
dicted by thin accretion disk models, as further dis- 
cussed in §4.5. 

4-1.1. The Far UV Continuum 

Zheng et al. (1996) have constructed a compos- 
ite quasar spectrum based on HST spectra of 101 
quasars at z > 0.33. They find a far UV (FUV) slope 
(1050A-350A) of (a FUV ) = -1.77±0.03 for RQQ and 
(o;fuv) = —2.16 ±0.03 for RLQ, with slopes of <~ — 1 
m the 2000A-1050A regime. The Zheng et al. mean 
spectra, presented in Fig. 6, together with the PSPC 
mean spectra, suggest that the FUV power-law con- 
tinuum extends to the soft X-ray band. In the case of 
RQQ there is remarkable agreement in both slope and 
normalization of the soft X-ray and the FUV power- 
law continua. RLQ are predicted to be weaker than 
RQQ at - 100 eV by both the FUV and the PSPC 
composites. It thus appears that there is no extreme 
UV sharp cutoff in quasars, and that the fraction of 
bolomctric luminosity in the FUV regime is signifi- 
cantly smaller than assumed. 

The UV to X-ray continuum suggested in Fig. 6 is 
very different from the one predicted by thin accre- 
tion disk models (§4.5), and suggested by photoion- 
ization models. In particular, it implies about four 
times weaker FUV ionizing continuum compared with 
the MF continuum which was deduced based on the 
He II A1640 recombination line. 

One should note, however, that the Zheng et al. 
sample is practically disjoint from our low z sample, 



so it may still be possible that low z quasars have a 
different FUV continuum. 

4.2. Intrinsic Absorption 

As shown in Fig. 2, the H I column deduced from 
our accurate 21 cm measurements is consistent for all 
objects with the best fit X-ray column. It is quite 
remarkable that even in our highest S/N spectra the 
21 cm and X-ray columns agree to a level of about 
1 x 10 19 cm~ 2 , or 5-7%. This agreement is remarkable 
since the 21 cm line and the PSPC are actually mea- 
suring the columns of different elements. Most of the 
soft X-ray absorption is due to He I or He II, rather 
than H I, and the H I column is indirectly inferred as- 
suming the column ratio H I/He 1= 10. The fact that 
the 21 cm line and the PSPC give the same H I column 
implies that the H I/He column ratio at high Galactic 
latitudes must indeed be close to 10. The dispersion 
in the H I/He column ratio is lower than 20% (based 
on typical quasars in our sample), and may even be 
lower than 5% (based on our highest S/N spectra). 
There is therefore no appreciable Galactic column at 
high Galactic latitudes where the ionized fraction of 
H differs significantly from the ionized fraction of He, 
as found for example in H II regions (e.g. Osterbrock 
1989). 

The consistency of the 21 cm and X-ray columns 
also indicates that the typical column of cold gas in- 
trinsic to the quasars in our sample must be smaller 
than the X-ray JVh uncertainty, or about 3 x 10 19 (1 + 
z) 3 cm~ 2 . An additional indication for a lack of an 
intrinsic cold column in quasars comes from the fact 
the the strong correlations of a x with the emission 
line parameters described above (§3.3) become weaker 
when we use a x from the free fit rather than a x 
from the fit with N§^. This indicates that A/gf is 
closer to the true Nn than the free -ZVh (see discus- 
sion in Paper I). In our highest S/N spectra we can 
set an upper limit of ~ 2 x 10 19 cm~ 2 on any intrin- 
sic absorption. As discussed in Paper I, the lack of 
intrinsic X-ray column for most quasars is consistent 
with more stringent upper limits set by the lack of a 
Lyman limit edge, as well as the He I and the He II 
bound-free edges in a few very high z quasars. 

Puchnarewicz et al. (1996) suggest that the strong 
a x vs. a ox correlation in their X-ray selected sam- 
ple is due to absorption of the optical and soft X-ray 
emission by cold gas and dust. They show that the a x 
vs. a ox correlation for the 10 objects in Paper I can 
be explained by a universal spectral shape absorbed 
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by a gas with a column of up to Nn— 3 x 10 20 cm~ 2 
(see their Figure 16). As described above, such ab- 
sorbing columns are clearly ruled out by our high S/N 
spectra. 

The mean a D in the Puchnarewicz et al. sample is 
—0.92 ± 0.07, which is significantly steeper than the 
mean a for optically selected quasars, e.g. a median 
of -0.2 for 105 PG quasars (Neugebauer et al. 1987), 
a median of —0.32 for 718 LBQS quasars (Francis et 
al. 1991), and (a Q ) = — 0.36 ± 0.05, in our sample. 
Puchnarewicz et al. suggested that the much flatter 
a Q of the PG quasars is a selection bias since these 
quasars were selected by the strength of their UV ex- 
cess. However, the PG sample was selected on the ba- 
sis of the color criterion U — B < —0.44, which using 
the flux transformations of Allen (1973), corresponds 
to a Q > —1.8. Thus, most of the red quasars dis- 
covered by Puchnarewicz et al. fit into the PG color 
criterion. The difference between the soft X-ray se- 
lected and optically selected quasars must reflect the 
true tendency of quasars selected above 0.5 keV to be 
significantly redder than optically selected quasars. 
These red quasars may very well be affected by a large 
absorbing column (-/V H > 10 21 cm -2 ), as suggested by 
Puchnarewicz et al. 

Intrinsic absorption is common in Seyfcrt 1 galax- 
ies. About half of the primarily X-ray selected Seyfert 
galaxies observed by Turner & Pounds (1989) using 
the EXOSAT LE+ME detectors, by Turner, George 
& Mushotzky (1993) using the ROSAT PSPC, and 
by Nandra & Pounds (1994) using the GINGA LAC 
for a largely overlapping sample, show low energy ab- 
sorption, or spectral features inconsistent with the 
simple power-law typically observed above 2 keV. 
Quasars are very different. Excess absorption pro- 
duces significant spectral features only in one object 
(PG 1114+445, see paper I), i.e. - 5% (1/23) of 
the objects, and the absorbing gas is partially ionized 
("warm"), rather than neutral. Given the typical S/N 
in our sample we estimate that a partially ionized ab- 
sorber which produces t > 0.3 can be ruled out in 
most of our objects. We cannot rule out partial ab- 
sorption, or complete absorption and scattering, by a 
very high column density (-/Vh> 10 24 cm~ 2 ) gas since 
such effects may only suppress the flux level with- 
out affecting the spectral shape, and without induc- 
ing significant spectral features. As described in §4.8, 
we suspect that such strong absorption may indeed 
be present in about <~ 10% (3/23) of the quasars in 
our sample (the X-ray weak quasars). 



4.3. Implications of the Continuum-Continuum 
Correlations 

The continuum-continuum luminosity correlations 
found here are all weaker than found in Paper I. This 
is mostly due to the three X-ray weak quasars which 
were not present in Paper I. For example, in Paper I 
we found that /0.3 kov can be predicted to within a 
factor of two, once /1.69 nm is given. This statement 
is still valid if the 4 extreme objects labeled in Fig. 3 
middle panel are excluded. The implications of the 
near IR versus X-ray luminosity correlation on the 
X-ray variability power spectrum were discussed in 
Paper I. 

In Paper I we noted the similarity ( 
— 1.50, which was also noted by Brunner et al. (1992) 
and Turner, George & Mushotzky (1993). However, 
we argued there that this similarity is only fortuitous, 
and that it does not imply that the X-ray power 
law can be extrapolated into the UV since the op- 
tical slope is significantly different. Here we find 
that the relation (a ox ) ~ (a x ) holds only roughly 
for the complete sample where (a ox ) = — 1.55 ± 0.24, 
and (a x ) = — 1.62 ± 0.45. This relation docs not 
hold when the sample is broken to the RQQ where 
(a ox ) = -1.56 ± 0.26, and (a x ) = -1.72 ± 0.41, 
and to the RLQ where (a ox ) — —1.51 ± 0.16, and 
(a x ) = -1.15 ±0.27. 

A significantly flatter (a ox ) is obtained when the 
three X-ray weak quasars, and the absorbed quasar 
PG 1114+445 are excluded. Thus, "normal" RQQ 
quasars in our sample have (a ox ) = — 1.48 ± 0.10, 
(a x ) = -1.69 ± 0.27, while for the RLQ (a ox ) = 
-1.44+0.12, (a x ) = -1.22 + 0.28, where the ± de- 
notes here and above the dispersion about the mean, 
rather than the error in the mean. 

The a x versus a ox correlation found here is weaker 
than in Paper I due to the presence of the X-ray weak 
quasars. However, the other 20 quasars appear to 
follow a trend of increasing a x with increasing a ox 
(Fig.5e), indicating as discussed in Paper I that a 
steep a x is generally associated with a weak hard X- 
ray component (at 2 keV), rather than a strong soft 
excess. The only object which clearly violates this 
trend is PG 1626+554 (Fig.3), which has both a steep 
a x and a strong soft excess. Puchnarewicz, Mason 
& Cordova (1994) and Puchnarewicz et al. (1995a; 
1995b) present PSPC spectra of three AGNs with ex- 
tremely strong soft excess, where Ln 2 krv > L 

to ! U.Z KOV 300Q A 

Our sample suggests that such objects are most likely 
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rare, as can also be inferred from the selection criteria 
of Puchnarewicz et al. who selected their three ob- 
jects from the ROSAT WFC all sky survey, in which 
only five AGNs were detected. This selection crite- 
rion implies that these AGNs must have a very high 
far UV flux. 

The soft X-ray selected quasars in the Puchnarewicz 
et al. (1996) sample have (a ox ) = -1.14 ± 0.02, 
and none of their quasars is "X-ray weak" , i.e with 
dox < —1.6. The absence of X-ray weak quasars in 
their sample is clearly a selection effect. The small 
survey area implies that most quasars in their sample 
are optically rather faint (rag ~ 18-19). "Normal" 
a ox quasars in their sample produce a few hundred 
PSPC counts, but "X-ray weak" quasars are below 
their detection threshold. The abundance of "X-ray 
loud" quasars (i.e a ox > —1) in the Puchnarewicz 
et al. sample is consistent with their rarity in op- 
tically selected samples. For example, quasars with 
ct x = — 1 are about 20 times fainter at 3000A than 
quasars with a ox = —1.5, for the same L x . Since the 
surface density of quasars increases as ~ /g ' (see 
§2.2.2.1 and Figure 1 in Hartwick & Schade), where 
/ B is the B band flux, there arc about 700 times more 
of these fainter quasars per B magnitude per square 
degree. Thus, even if only 0.3% of quasars at a given 
/b have a ox = — 1, there would still be twice as many 
quasars with a ox = — 1 than a ox = — 1.5 per square 
degree in an X-ray flux limited sample, such as the 
Puchnarewicz et al. sample. 

4.4. Implications of The Continuum-Line Cor- 
relations 

The presence of the strong correlations of a x with 
the H/3 FWHM, with L [Q nl] , and with the Fe II/H/3 
ratio described in Paper I is verified here. The corre- 
lation coefficients for the complete sample are compa- 
rable or somewhat smaller than those found in Paper 
I, but since the sample is larger the significance level 
is now much higher (Fig. 5). We also report here an 
additional strong correlation of a x with the ratio of 
[O III] peak flux to the H/3 peak flux, which is one 
of the emission line parameters measured by Boro- 
son & Green. The a x versus H/3 FWHM correlations 
is the strongest correlation we find between any of 
the X-ray continuum emission parameters and any of 
the emission line parameters reported by Boroson & 
Green (with the addition of line luminosities reported 
in Table 5). 

The a x versus H/3 FWHM correlation is much 



stronger than the well known a x correlation with ra- 
dio loudness (rg = 0.26, Pr=0.23 in our sample, but 
see Wilkes & Elvis 1987 and Shastri et al. 1993 for 
stronger correlations). Thus, the fact that the average 
a x in RLQ is significantly flatter than in RQQ (§4.1) 
may be completely unrelated to the presence of radio 
emission, it may just reflect the fact that RLQ tend 
to have broader lines than RQQ (e.g. Tables 3 and 

5 in Boroson & Green). This appears to be the case 
in our sample, where the RLQ follow the same a x 
versus H/3 FWHM distribution defined by the RQQ 
(Fig. 5a). This intriguing suggestion can be clearly 
tested by comparing a x for RLQ and RQQ of similar 
H/3 FWHM. 

Boiler, Brandt & Fink (1996) studied in detail nar- 
row line Seyfert 1 galaxies (NLS1) and they also find 
an apparently significant trend of increasing a x with 
increasing H/3 FWHM. However, the scatter in their 
sample is significantly larger than that found here. 
In particular they find a large range of a x for H/3 
FWHM< 2000 km s _1 , where only a few objects are 
available in our sample. The overall larger scatter in 
the Boiler, Brandt & Fink data is probably due in 
part to the generally larger statistical errors in their 
a x determinations. Large systematic errors may also 
be induced by the use of H/3 FWHM from a variety of 
sources. The measured H/3 FWHM can be sensitive 
to the exact measuring procedure, such as continuum 
placement, subtraction of Fe II blends, and subtrac- 
tion of the narrow component of the line (produced 
in the narrow line region) which may not be well re- 
solved in low resolution spectra. For example, Shastri 
et al. 1993 and Boroson & Green measured the H/3 
FWHM independently for 13 overlapping objects, in 8 
of which their values differ by more than 1000 km s _1 . 
Other than these technical reasons the increased scat- 
ter may represent a real drop in the strength of the 
correlation when the luminosity decreases from the 
quasar level studied here to the Seyfert level studied 
by Boiler, Brandt & Fink. One should also note that 
intrinsic absorption is probably common in Seyfert 
1 galaxies (Turner & Pounds, 1989; Turner, George 

6 Mushotzky, 1993), and such an absorption may 
lead to a large systematic error in a x unless a high 
S/N spectrum is available indicating that features are 
present. 

Wang, Brinkmann & Bergeron (1996) analyzed 
PSPC spectra of 86 AGNs, including 22 of the 23 
quasars from our sample. Their sample is more 
heterogeneous than ours and includes some high z 
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quasars and a number of AGNs selected by their 
strong Fe II emission. The various correlations found 
by Wang et al. are typically similar, or somewhat 
weaker than found here. For example, their(our) val- 
ues are r s = -0.73(-0.79) for a x versus H/3 FWHM, 
and r s = 0.65(0.714) for a x versus Fe II/H/3. The 
somewhat smaller values found by Wang et al. may 
result from their inclusion of z > 0.4 quasars, where 
a x measures a higher energy slope than measured 
here, and from their use of a free Nh fit (limited from 
below by iVjfj ), which increases the random error in 
a x (see Paper I). 

We verify the strong correlation between Ljj/3 and 
L2 keV found in Paper I. The correlations of the other 
lines with X-ray luminosity are significantly weaker 
than found in Paper I, and they are only marginally 
significant. Corbin (1993) found significant correla- 
tions of L 2 kov with Fe II/H/3 (r s = -0.474), and 
of L2 kcv with the H/3 asymmetry (rs = —0.471). 
We find that neither correlation is significant in our 
sample (r s = -0.288, Pr=0.19, and r s = -0.106, 
Pr=0.63). Since we can only test for correlation with 
rs > 0.64, we cannot securely exclude the presence of 
the correlations reported by Corbin. 

As discussed in Paper I, the a x versus L[ G in] cor- 
relation can be used to place a limit on the a x vari- 
ability on timescales shorter than a few years. Given 
the scatter in this correlation we estimate that a x 
should not vary by significantly more than 0.3 on 
these timescales. 

It is hard to interpret the a x versus [O III] to H/3 
peak flux ratio correlation since the physical mean- 
ing of the [O III] to H/3 peak flux ratio parameter 
defined by Boroson & Green is rather obscure. The 
[O III] peak flux is related to the width of [O III], 
and the [O III] to H/3 peak flux ratio may thus partly 
reflect the FWHM ratio of these lines. Thus, this 
correlation may represent a correlation of the [O III] 
FWHM with a x . High spectral resolution measure- 
ments of the [O III] line profile are required to test 
this possibility. 

4.5. Inconsistency with Thin Accretion Disk 
Models 

Figure 7 presents the continuum emission from two 
thin accretion disk models. The models are for a disk 
around a rotating black hole, and viscous stress which 
scales like the ^/P gas P rac i, where P ra( j is the radiation 
pressure and P gas is the gas pressure (Laor & Netzer 



1989). Significant soft X-ray emission is obtained for 
disks with a high accretion rate and a small inclina- 
tion. However, as discussed by Fiore et al. (1995), 
the observed soft X-ray spectral slope is always much 
flatter than the one produced by a thin 'bare' accre- 
tion disk model. As noted above there is no indication 
in the 0.2-2 keV band for a very steep and soft "ac- 
cretion disk" component. 

Although thin disks cannot reproduce the 0.2- 
2 keV spectral shape, they may still be able to con- 
tribute a significant fraction of the flux at the lowest 
observed energy, i.e. 0.2-0.3 keV, above which a non 
thermal power-law component sets in. As noted by 
Walter et al. (1994) and in Paper I, accretion disk 
models predict a large dispersion in the optical/soft 
X-ray flux ratio, and the strong correlation between 
these fluxes argues against the idea that a thin disk 
produces both the optical and soft X-ray emission. 
The arguments put by Walter et al. and in Paper I 
were only qualitative, and were not based on actual 
disk models. Furthermore, the objects in the small 
sample of Walter et al. were selected from known op- 
tically bright AGNs, and they also had to be bright 
soft X-ray sources since most spectra were obtained 
from the ROSAT all sky survey. Thus, these objects 
were a priori selected to be bright at both optical- UV 
and at soft X-rays, and the absence of a large scat- 
ter in the UV/soft X-ray flux ratio may just reflect 
the sample selection criteria. Such selection effects 
are not present in our sample since the sample was 
defined independently of the X-ray properties, and 
X-ray spectra were obtained for all objects. 

Below we describe a detailed calculation of the ex- 
pected distribution of optical/soft X-ray flux ratio 
for a complete optically selected sample based on the 
thin disk models of Laor & Netzer (1989), and show 
that such models cannot be reconciled with the ob- 
served distribution of optical/soft X-ray flux ratio in 
our complete sample. 

The optical/soft X-ray flux ratio, a OSl of a given 
disk model depends on the black hole mass, accretion 
rate m, and inclination angle 9 = cos^ 1 fi. We now 
need to determine what distributions of these param- 
eters will be consistent with the observed luminosity 
function in a complete optically selected sample. The 
intrinsic distribution of disk inclinations must be ran- 
dom. However, the observed distribution depends on 
the shape of the luminosity function of quasars, and 
possible obscuration effects, as described below. 

The luminosity function of quasars is parametrized 
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using the number density of quasars per unit volume 
per magnitude $ = d 2 N/dMdV, and it is well fit by 
a power-law over a restricted range of magnitude, M. 
Using Figure 2 in Hartwick & Schade (1990) we find 
log $ = 0.55M + c for z < 0.2 and log $ = 0.66M + c 
for 0.4 < z < 0.7, where c is a constant. Since our 
sample is restricted to z < 0.4 we assume log$ = 

0. 6M + c. Using the relation M = — 2.5 log L + c we 
get that dn/dL cx L~ 2 ' 5 , where n = dN/dV. 

The apparent luminosity of a flat disk L app is re- 
lated to its intrinsic luminosity through L app = 2uL, 
neglecting limb darkening effects which steepen the 
/i dependence, and relativistic effects which flatten 
the \i dependence. This provides a reasonable ap- 
proximation in the optical-near UV regime (see Laor, 
Nctzer, & Piran 1990). Assuming dN/d/j, = const, 

1. e. a uniform distribution of inclination angles for 
the intrinsic quasar population, we would like to find 
dn/du for a given L app . When L app is fixed, /i cx L , 
and substituting u in the expression for dn/dL we get 
dn/du oc u - 5 . Thus although the disks are assumed 
to have a uniform distribution of inclination angles 
the observed distribution at a given L app is biased 
towards face on disks. 

To reproduce the observed luminosity function we 
choose two values of to = 0.1,0.3, where to is mea- 
sured in units of the Eddington accretion rate. Since 
L a m 9 m, where mg is the black hole mass in units of 
10 9 Mq the required mass distribution is dn/dmg oc 
nig 2 ' 5 . The observed number of objects in a flux 
limited sample is dN ^/dL oc dn/dL x V(L), where 
V{L) cx L 3 / 2 is the observable volume for a flux lim- 
ited sample, such as the BQS sample. We therefore 
select a mass distribution of dN b/dnig cx nig 1 . 

Figure 8 compares the observed distribution of a os , 
as a function of vL v at 3000A, with the one expected 
from thin accretion disk models with the parameter 
distribution described above. Thin disk models can- 
not account for the very small scatter in ct os . 

The range of observed disk inclinations may ac- 
tually be smaller than assumed here. For example, 
for a certain range of inclinations the disk may be 
completely obscured by an optically thick torus, as 
suggested in unification schemes for RQQ (e.g. An- 
tonucci 1993). However, even if u is fixed at a given 
value for all AGNs (say u — 1 which corresponds to 
the points extending from a os = —1.5 on the left axis 
to \ogvL v = 46.5 on the bottom axis of Fig. 8), the 
range in to 9 and to will stil produce a range in a os 
which is much larger than observed. 



The X-ray power-law emission is most likely pro- 
duced by Comptonization of the thermal disk emis- 
sion in a hot corona above the disk (e.g. Czerny 
& Elvis 1987). The slope and normalization of the 
power-law component are determined by the tem- 
perature and electron scattering optical depth in the 
corona (e.g. Sunyaev & Titarchuk 1985; Titarchuk 
& Lyubarskij 1995). The small range in a os implies 
that some physical mechanism which couples the op- 
tical and soft X-ray emission processes must be oper- 
ating, e.g. through a feedback which regulates both 
the temperature (see Haardt & Maraschi 1993) and 
the optical depth of the corona. 

As pointed out by various authors (Ross, Fabian & 
Mineshige, 1992; Shimura & Takahara 1995; Dorrer 
et al. 1996), and shown in Fig. 7, simple thin accretion 
disks with no corona can produce a significant flux be- 
low 1 keV. For various disk model parameters a os can 
in fact be significantly flatter than observed (Fig. 8), 
yet such extreme flat optical-soft X-ray spectra are 
only rarely observed (e.g. Puchnarewicz 1995a). The 
flattest spectra are expected for disks which are close 
to edge on (e.g. Laor, Netzer & Piran 1990), and one 
therefore needs to assume that such disks are not ob- 
servable. This is indeed expected in AGNs unification 
schemes which invoke obscuring material close to the 
disk plane. Alternatively, the accreted material may 
form a geometrically thick, rather than a thin, con- 
figuration close to the center, which would display a 
smaller inclination dependence. 

4.6. Inconsistency with Optically Thin Free- 
Free Emission Models 

As was clearly demonstrated by Fiore et al. (1995) 
for 6 low redshift quasars with a high S/N PSPC spec- 
tra, and by Sicmiginowska et al. (1995) using EIN- 
STEIN data for 47 quasars from Elvis et al. (1994), 
isothermal optically thin pure free-free emission mod- 
els (Barvainis 1993) cannot fit the observed UV to soft 
X-ray energy distribution in AGNs. Furthermore, as 
was pointed out by Kriss (1995), and Hamman et al. 
(1995), optically thin free-free emission can also be 
ruled out based on the observed UV line emission. 

4.7. On the origin of the a x versus H/3 FWHM 
correlation 

What is the physical process behind the a x versus 
H/3 FWHM correlation? In Paper I we speculated 
that this may either be an inclination effect, or that 
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it could be an L/LEdd effect. Various authors raised 
the interesting suggestion that steep a x quasars may 
be analogous to 'high'-state Galactic black hole can- 
didates (e.g. White, Fabian & Mushotzky 1984; Fiore 
& Elvis 1995; Pounds, Done & Osborne 1995), which 
display a steep slope in the soft and the hard X-ray 
bands when their brightness increases. The physi- 
cal interpretation for this effect described by Pounds, 
Done & Osborne (1995) is that the hard X-ray power- 
law is produced by Comptonization in a hot corona 
and that as the object becomes brighter in the optical- 
UV, Compton cooling of the corona increases, the 
corona becomes colder, thus producing a steeper X- 
ray power-law. This is obviously far from being a 
predictive model since the corona heating mechanism 
is not specified, and it is implicitly assumed that 
the coronal heating does not increase much as the 
quasar becomes brighter. However, the narrow H/3 
line profiles provide independent evidence that steep 
a x quasars may indeed have a higher L/LEdd, as fur- 
ther described below. 

The L/LEdd of quasars can be estimated under two 
assumption: 1. The bulk motion of the gas in the 
broad line region is virialized, i.e. Av ~ yj GM/r, 
where Av=H/3 FWHM. This gives 

Awaooo = 2.19ml /2 Ro.\ /2 , 

where Aw 300 o = Aw/3000 km s _1 , m 9 = M/10 9 M , 
and Ro.i — R/0.1 pc. 2. The size of the broad 
line region is determined uniquely by the luminosity, 
Ro.i = L^ 2 pc, where L 46 = L Bo iomotric/10 46 . This 
scaling is consistent with reverberation line mapping 
of AGNs (Peterson 1993; Maoz 1995), and is theoret- 
ically expected if the gas in quasars is dusty (Laor 
& Draine 1993, Netzer & Laor 1994). Combining as- 
sumptions 1 and 2 gives 

Av 300 o = 2.19mg /2 L 46 1/4 , 
and thus the mass of the central black hole is 

m 9 = 0.21Au| 000 L^ 2 . 
Using LEdd,46 = 12.5mg one gets 

L/LEdd = 0.38A W30 2 L^ 2 . 

Thus, given the two assumptions made above, narrow 
line quasars should indeed have a high L/LEdd, as 
previously suggested based only on their steep a x , 
and analogy to Galactic black hole candidates. 



To test whether L/Lsdd is indeed the underlying 
parameter which determines a x , rather than just the 
H/3 FWHM, we looked at the correlation of a x ver- 
sus AwjqqqL^ 2 displayed in Fig.5f, where we used 
the 1.7 //m luminosity and the relation LBoiomotric = 
I5L14.25 (see Fig. 7 in Laor & Draine). This corre- 
lation is not as strong as the a x versus H/3 FWHM 
correlation, but it certainly appears suggestive. Note 
that L/L E dd > 1 for some of the objects in the sam- 
ple. These values are well above the thin accretion 
disk limit (L/L^dd = 0.3, Laor & Netzer 1989) and 
suggest a thick disk configuration. However, the as- 
sumptions used above to infer L/LEdd are more quali- 
tative than quantitative since both the luminosity and 
the velocity field in the broad line region may not be 
isotropic and therefore the presence of L/LEdd > 1 
cannot be securely deduced. 

The Pounds et al. mechanism implies that a steep 
a x is associated with a weak hard X-ray component, 
and as described in §3.2, this indeed appears to be the 
trend in our sample. If the Pounds et al. mechanism 
is true then steep a x AGNs should also have a steep 
hard X-ray power-law. We are currently pursuing this 
line of research using ASCA and SAX observations of 
our sample. 

An additional hint that a steep a x may indeed be 
associated with a high L / L^dd comes from the anec- 
dotal evidence described by Brandt, Pounds & Fink 
(1995), Brandt et al. (1995), Grupc et al. (1995), and 
Forster & Halpern (1996) where a number of Seyfert 
galaxies with a steep a x display rapid, large ampli- 
tude, soft X-ray variability, which as Boiler, Brandt, 
& Fink discuss may imply a low mass black hole, and 
thus a high L/L Edd . We are currently pursuing a 
more systematic study of the soft X-ray variability 
properties of broad versus narrow line quasars using 
the ROSATHRI. 

The Pounds et al. suggestion is very appealing 
since it allows a physical explanation for the tight 
correlation of apparently completely unrelated quan- 
tities. Although it is not clear a priori that a x must 
steepen with increasing L/LEdd, it appears that this 
is indeed what happens in Galactic black hole candi- 
dates. 

Wang et al. also suggested that steep a x objects 
have a high L/L^dd based on the fact that the frac- 
tion of luminosity emitted in the X-ray regime in thin 
accretion disk models increases with L/LEdd, as dis- 
cussed above in §4.5. However, if this were indeed the 
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physical process behind the a x versus H/3 FWHM cor- 
relation then one would expect high L/L^dd objects 
to have a high soft X-ray to optical flux ratio, while 
we find no correlation between H/3 FWHM and a os 
(r s = -0.079, Table 5). 

We note in passing that one does not need to elim- 
inate the normal broad line region in narrow line 
AGNs, as one of the options suggested by Boiler et 
al., and Pounds et al. The lines are narrow simply 
because of the lower black hole mass. The broad line 
emitting gas does not extend much closer to the cen- 
ter in narrow line AGNs, as it does not extend much 
closer to the center in other AGNs, simply because 
of the effects of a higher ionization parameter, and a 
higher gas density, each of which quenches line emis- 
sion. 

4.8. The X-ray Weak Quasars 

Two of the quasars in our sample, the RQQ PG 1001 
+054 and PG 1411+442, and possibly also the RLQ 
PG 1425+267 appear to form a distinct group which 
we term here "X-ray weak" quasars, where the nor- 
malized X-ray luminosity is a factor of 10-30 smaller 
than the sample median. The position of these 
quasars as outliers can be noticed in the near IR nor- 
malized flux distribution (Fig. 3), in the a x versus a ox 
correlation (Fig.5e), and in the H/3 versus 2 keV and 
0.3 keV luminosity correlations (Fig. 4). The first two 
indicators are based on the spectral shape, but the 
last one is independent of the spectral shape, and it 
also suggests a deficiency of the X-ray luminosity by 
a factor of 10-30 relative to the one expected based 
on the H/3 luminosity. An apparently bimodal distri- 
bution in a ox can also be seen in Figure 5b of Wang 
et al. where 6 of their 86 quasars appear to form a 
distinct group with a ox < — 2. 

No bimodality of a ox is seen in the Avni & Tanan- 
baum (1986) EINSTEIN study of the PG quasars. 
All the "X-ray weak" quasars found by the ROSAT 
PSPC have a ox < —2, but none of the quasars de- 
tectd by Tananbaum et al. have a ox < —2 (see Fig. 8 
in Avni & Tananbaum). The lack of a ox < — 2 and 
bimodality in the Tananbaum et al. sample probably 
reflects its incompleteness, as only 86% of the quasars 
they observed were detected. 

Although the three X-ray weak quasars in our sam- 
ple stand out in luminosity correlations, they con- 
form well to the a x correlations (Figs.5a-d). They 
thus have the "right" slope but the "wrong" flux 



level. Why are these quasars different? A simple 
answer is that for some unknown reason the X-ray 
emission mechanism, most likely Comptonization by 
T > 10 8 K electrons, tends to be bimodal, and in 
about 10% of quasars (or in all quasars for <~ 10% of 
the time) the X-ray flux level is strongly suppressed, 
while the spectral slope is not affected. Another op- 
tion is that these are just normal quasars where the 
direct X-ray flux happens to be obscured. In this case 
what we see is only the scattered X-ray flux. Pho- 
toionization calculations indicate that a few percent 
of the direct flux will be scattered, depending on the 
covering factor of the absorber and the ionization pa- 
rameter. If the ionization parameter is large enough 
then the scattering will be mostly by free electrons 
which preserves the spectral shape (see Netzer 1993, 
and by Krolik & Kriss 1995). Such scattering will ex- 
plain why the flux level is strongly reduced, while the 
spectral shape is not affected. Note that the obscur- 
ing matter should be transparent in the visible range, 
as is the case with the absorbing matter in BALQSO. 

Additional hints towards this interpretation come 
from the fact that PG 1411+442 is a broad absorp- 
tion line quasar (BALQSO, Malkan Green & Hutch- 
ings, 1987), and the UV absorbing gas may also pro- 
duce soft X-ray absorption, as may also be the case in 
PHL 5200 (Mathur Elvis & Singh, 1995). In addition, 
Green & Mathur (1996) find that BALQSO observed 
by R OS A T have a ox <^ —1.8, i.e. as observed here 
in the X-ray weak RQQ. Another hint is provided 
by the fact that PG 1114+445 is also somewhat un- 
derluminous at 0.3 keV (Fig. 3b), and this quasar is 
most likely seen through a warm absorber. The X-ray 
weak quasars could therefore be more extreme cases 
of PG 1114+445 and have an absorbing column which 
is large enough to completely absorb the direct soft 
X-ray emission. 

Note that PG 1425+267 is a RLQ, while all BALQSO 
are known to be RQQ (Stocke et al. 1992). It would 
thus seem implausible to suggest that PG 1425+267 is 
a BAL. However, PG 1425+267 has about the same 
relatively steep a ox , compared to other RLQ, as in 
3C 351 (Fiore et al. 1993), where X-ray absorption 
by warm gas is observed together with resonance UV 
absorption lines (Mathur et al. 1994) which are nar- 
rower than in 'proper' BALQSO. 

Forthcoming HST spectra of all 23 quasars in our 
sample will allow us to test if there is a one to one 
correspondence between X-ray weakness and broad 
absorption lines, i.e. if all X-ray weak quasars are 
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BALQSO, and not just that all BALQSO are X-ray 
weak, as strongly suggested by the Green & Mathur, 
and the Green et al. (1996) results. 

A simple test of whether these are truly "X-ray 
weak quasars" , or just normal highly absorbed quasars, 
can be done by looking at their hard X-ray emis- 
sion. If the X-ray column is below 10 24 cm~ 2 then 
the obscuring material would become transparent at 
E < 10 keV, and the observed hard X-ray emission 
will rise steeply above the cutoff energy, as seen in var- 
ious highly absorbed AGNs, such as Mkn 3 (Iwasawa 
ct al. 1994), NGC 5506 (Nandra & Pounds 1994), 
NGC 6552 (Fukazawa, et al. 1994; Reynolds et al. 
1994), and NGC 7582 (Schachter et al. 1996). One 
may also expect significant spectral features produced 
by the obscuring material (e.g. Matt et al. 1996), 
depending on the ionization state of this material. 
Forthcoming ASCA observations of PG 1411+442 
and PG 1425+267 will allow us to test this scenario. 

Another prediction is that the X-ray weak quasars 
should show lower variability compared with other 
quasars of similar X-ray luminosity. This is because: 
1). they are intrinsically more X-ray luminous, and 
variability amplitude tends to drop with increasing 
luminosity (Barr & Mushotzky 1986; also Fig. 9 in 
Boiler, Brandt & Fink). 2). the scattering medium 
must be significantly larger than the X-ray source, 
and short time scale variability will be averaged out. 
If the X-ray weak quasars are just due to large am- 
plitude intrinsic variability of the soft X-ray emission, 
as seen in some steep narrow line Seyfert 1 galaxies 
(§4.7), then one may expect the exact opposite be- 
havior, i.e. these quasars may become significantly 
brighter at soft X-rays at some stage in the future. 

5. SUMMARY 

We defined a complete sample of 23 optically se- 
lected quasars which includes all the PG quasars at 
z < 0.400, and 7Vg a /< 1.9 x 10 20 cm" 2 . Pointed 
-RO^TPSPC observations were made for all quasars, 
yielding high S/N spectra for most objects. The high 
quality of the ROSAT spectra allows one to determine 
the best fitting a x with about an order of magnitude 
higher precision compared with previously available 
X-ray spectra. In this paper we report the observa- 
tions of 13 quasars not described in Paper I, analyze 
the correlation of the X-ray properties of the complete 
sample with other emission properties, determine the 
mean X-ray spectra of low z quasars, discuss the pos- 



sible origin of the a x versus H/3 FWHM correlation, 
the nature of X-ray weak quasars, and the physical 
origin of the soft X-ray emission. Our major results 
are the following: 

1. The spectra of 22 of the 23 quasars are consis- 
tent, to within ~ 10 — 30%, with a single power- 
law model over the rest frame range 0.2 — 2 keV. 
There is no evidence for significant soft excess 
emission with respect to the best fit power-law. 
We place a limit of <~ 5 x 10 19 cm~ 2 on the 
amount of excess foreground absorption by cold 
gas in most of our quasars. The limits are 
<~ 1 x 10 19 cm~ 2 in the two highest S/N spectra. 

2. Significant X-ray absorption by partially ionized 
gas ( "warm absorber" ) in quasars is rather rare, 
occurring for <; 5% of the population, which 
is in sharp contrast to lower luminosity AGNs, 
where significant absorption probably occurs for 
<~ 50% of the population. 

3. The average soft X-ray spectral slope for RQQ 
is (a x ) = —1.72 + 0.09, and it agrees remarkably 
well with an extrapolation of the mean 1050A- 
350A continuum recently determined by Zheng 
ct al. (1996) for z > 0.33 quasars. For RLQ 
(a x ) = —1.15 ± 0.16, which suggests that RLQ 
quasars are weaker than RQQ below 0.2 keV, as 
suggested also by the Zheng et al. mean RLQ 
continuum. These results suggest that there is 
no steep soft component below 0.2 keV. 

4. Extensive correlation analysis of the X-ray con- 
tinuum emission parameters with optical emis- 
sion line parameters indicates that the strongest 
correlation is between a x , and the H/3 FWHM. 
A possible explanation for this remarkably strong 
correlation is a dependence of a x on L/L^dd, as 
observed in Galactic black hole candidates. 

5. There appears to be a distinct class of "X-ray 
weak" quasars, which form ~ 10% of the pop- 
ulation, where the X-ray emission is smaller by 
a factor of 10-30 than expected based on their 
luminosity at other bands, and on their H/3 lu- 
minosity. 

6. Thin accretion disk models cannot reproduce 
the observed 0.2-2 keV spectral shape, and they 
also cannot reproduce the tight correlation be- 
tween the optical and soft X-ray emission. 
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7. The H I/He I ratio in the ISM at high Galac- 
tic latitudes must be within 20%, and possibly 
within 5%, of the total H/He ratio. 

The main questions raised by this study are: 

1. What is the true nature of X-ray quiet quasars? 
Are these quasars indeed intrinsically X-ray 
weak, or are they just highly absorbed but oth- 
erwise normal quasars? 

2. What physical mechanism is maintaining the 
strong correlation between the optical-UV and 
the soft X-ray continuum emission, or equiva- 
lcntly, maintaining a very small dispersion in 
the maximum possible far UV cutoff tempera- 
ture? 

3. What is the physical origin for the strong corre- 
lations between a x , and L[ G in], Fe II/H/3, and 
the peak [O III] to H/3 flux ratio? 

4. Is the soft X-ray emission indeed related to the 
presence of radio emission, or is it just a spuri- 
ous relation and the primary effect is related to 
the H/3 line width? Or, put differently, do RLQ 
and RQQ of similar H/3 FWHM have similar 

Extensions of the ROSAT PSPC survey described 
in this paper to the hard X-ray regime with ASCA 
and SAX, to the UV with HST, and soft X-ray vari- 
ability monitoring with the ROSAT HRI, which are 
currently being carried out, may provide answers to 
some of the questions raised above. These studies 
will also allow us to: 1) Test if steep a x quasars 
have a steep 2-10 keV slope, as expected based on 
the Pounds et al. L/L^dd interpretations. 2) Test if 
soft X-ray variability is indeed strongly tied to the H/3 
FWHM, as expected if the H/3 FWHM is an indica- 
tor of L/LEdd- 3) Explore the relation of the UV line 
emission properties to the ionizing spectral shape. 
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